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Abstract. We show that phase transitions in Ising systems with planar 
defects, i.e., disorder perfectly correlated in two dimensions are destroyed by 
smearing. This is caused by effects similar to but stronger than the Griffiths 
phenomena: Exponentially rare spatial regions can develop true static long- 
range order even when the bulk system is still in its disordered phase. Close 
to the smeared transition, the order parameter is very inhomogeneous in space, 
with the thermodynamic (average) order parameter depending exponentially on 
temperature. We determine the behavior using extremal statistics, and we 
illustrate the results by computer simulations. 



1. Introduction 

Phase transitions in random systems remain one of the important and only partially 
solved problems in statistical physics. Originally, it was thought that even weak 
disorder destroys any critical point because the disordered system divides itself up 
into spatial regions which undergo the transition at different temperatures, leading to 
a smeared transition in which there would be no sharp singularities in thermodynamic 
quantities (see pQ and references therein). However, it was soon found that classical 
continuous phase transitions are generically sharp in the presence of weak, short-range 
correlated disorder. If a clean critical fixed point (FP) fulfills the Harris criterion [2] 
v >2 /d, where v is the correlation length exponent and d is the spatial dimensionality, 
the disorder decreases under coarse graining. The system becomes asymptotically 
homogeneous at large length scales. Thus, the clean critical point is perturbatively 
stable against disorder, and the critical behavior of the random system is identical to 
that of the corresponding clean system. Macroscopic observables are self-averaging at 
the critical point, i.e., the relative width of their probability distributions goes to zero 
in the thermodynamic limit 00]. 

Even if the Harris criterion is violated, the transition will generically remain sharp, 
but the critical behavior will be different from that of the clean system. Depending 
on the fate of the disorder under coarse graining it can either be of finite or of infinite 
disorder type. In the first case the system remains inhomogeneous at all length scales 
with the relative strength of the inhomogeneities approaching a finite value for large 
length scales. These transitions are controlled by renormalization group (RG) fixed 
points with finite disorder. Macroscopic observables are not self-averaging, the relative 
width of their probability distributions approaches a size- independent constant 00]. 
Examples of critical points in this class include the dilute three-dimensional Ising 
model or classical spin glasses. The second case occurs when the relative magnitude of 
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the inhomogeneities increases without limit under coarse graining. The corresponding 
fixed points are called infinite-disorder fixed points. Here, the probability distributions 
of macroscopic variables become very broad (on a logarithmic scale) with the width 
increasing with system size. This can lead to activated rather than the usual power- 
law scaling at the critical point. A famous classical example of an infinite-disorder 
critical point occurs in the McCoy- Wu model 0. 

In recent years, another aspect of phase transitions in disordered systems, the 
Griffiths phenomena, has regained a lot of attention. Griffiths phenomena are non- 
perturbative effects of rare disorder fluctuations in the vicinity of a phase transition. 
They can be understood as follows: Disorder in general decreases the critical 
temperature T c from its clean value Tj?. In the temperature region T c < T < T c ° 
the system does not display global order. However, in an infinite system one will 
find arbitrarily large regions that are devoid of impurities with a small but non- 
zero probability. These rare regions or 'Griffiths islands' are locally in the ordered 
phase while the bulk system is still in the disordered phase. Their fluctuations are 
very slow because flipping them requires changing the order parameter in a large 
volume. Griffiths [S] was the first to show that this leads to a non-analytic free 
energy everywhere in the region T c < T < T®, which is now known as the Griffiths 
phase 7 or the Griffiths region. In generic classical systems, the contribution of 
the Griffiths singularities to thermodynamic (equilibrium) observables is very weak 
since the singularity in the free energy is only an essential one (HJ|S]. In contrast, the 
consequences for the dynamics are much more severe with the rare regions dominating 
the behavior for long times. In the Griffiths region, the spin autocorrelation function 
C(t) decays very slowly with time t, like lnC(t) ~ — (lnt) d > ( d-1 J for Ising systems 
|7IB?lllUll 11112) . and like lnC(t) ~ — t 1 / 2 for Heisenberg systems |1 11113] . More recently, 
these results have been confirmed by rigorous calculations for the equilibrium |14U15| 
and dynamic |16| properties of disordered Ising models. 

In many real systems, the disorder is not generated by point-like defects but by 
dislocations or disordered layers or grain boundaries. Often, extended impurities in 
a ci-dimensional system can be modeled by disorder perfectly correlated in dc = 1 or 
2 dimensions but uncorrelated in the remaining d± = d — dc dimensions. While it 
is generally accepted that long-range disorder correlations increase the effects of the 
impurities, their influence on the phase transition has been controversial. Early RG 
work E| based on a single expansion in e = 4 — d did not produce a critical fixed point. 
This was interpreted as a smeared transition or a first order one However, 
an additional expansion in the number dc of correlated dimensions |2UII21ll^] cured 
this problem and gave rise to critical fixed points with conventional power-law scaling. 
Further support for the sharp transition scenario came from Monte-Carlo simulations 
of a 3D Ising model with planar defects J22|. All the perturbative RG studies miss, 
however, the effects of rare regions and Griffiths phenomena. These effects have been 
studied in detail in the already mentioned McCoy- Wu model |S] , a disordered 2D Ising 
model in which the disorder is perfectly correlated in one dimension and uncorrelated 
in the other. The transition in this model was originally thought to be smeared. It was 
later found that it is sharp but of infinite-disorder type |24II25) . Similar behavior has 
been found in several random quantum systems which are related to classical systems 
with linear defects 1213122112111211112312] • Based on these results it was believed for 
a long time that phase transitions generically remain sharp even in the presence of 
extended disorder. 

In this paper we show, however, that this belief is not generally true. Specifically, 
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we show that for Ising order parameter symmetry, planar defects destroy the sharp 
continuous phase transition via effects which are similar to but stronger than the 
usual Griffiths phenomena: In systems with planar defects true static long-range 
order can develop on isolated rare regions. As a result, the order parameter is very 
inhomogeneous in space close to the smeared transition, and the thermodynamic 
order parameter depends exponentially on temperature. We also show that this 
disorder-induced smearing generically occurs in an entire class of transitions in systems 
with extended disorder. The paper is organized as follows. In Sec. 13 the model is 
introduced, and the mechanism for the smearing of the transition is explained. In Sec. 
121 we use Lifshitz tail arguments to determine the behavior in the 'tail' of the rounded 
transition, i.e., for a very low concentration of ordered islands. Section^is devoted to 
a numerical study of a three-dimensional infinite-range Ising model which illustrates 
the smearing of the phase transition. Finally, in Sec. we discuss the generality of the 
disorder-induced smearing, the relation to Griffiths phenomena, favorable conditions 
for observing the smearing, and the influence of the order parameter symmetry. 

2. Rare regions, inhomogeneous order, and smearing 

For definitcness we consider a d-dimensional 4 theory with random- T c type disorder 
completely correlated in the dc — 2 directions x\ and X2 but uncorrelated in the 
remaining dj_ = d — dc directions 2:3, . . . , x^. The action is given by 




Here 0(r) is a scalar field, t is the dimensionless distance from the clean critical 
point and 6t(r±) introduces the quenched disorder, rj_ is the projection of r on the 
uncorrelated directions Z3, ... ,374. We consider two different types of disorder. In the 
first type, Poisson (or dilution) disorder, 5t(r±) is given by 

St(r ± )=J2V[r ± -r ± (i)] (2) 

i 

where rj_(i) are the random positions of planar impurities of spatial density c, and 
V^(rj_) is a non-negative short-ranged impurity potential. The second type of disorder 
is a Gaussian distribution with zero mean and a correlation function 

(5t(r 1 _)St(v ± ')) = A 2 S(r ± -r s _') . (3) 

Now consider the effects of rare disorder fluctuations in this system: Analogous to 
the Griffiths phenomena, there is a small but finite probability for finding large spatial 
regions in rj_ direction which are more strongly coupled than the bulk system. For 
Poisson disorder, these are regions devoid of impurities while for Gaussian disorder it 
would be regions with a negative average (St(r±)). These rare regions can be locally in 
the ordered phase, even if the bulk system is still in the disordered phase. For Poisson 
disorder which fulfills 6t(r±) > this starts to happen right below the clean transition, 
i.e., for t < j22 while for the unbounded Gaussian disorder, locally ordered rare 
regions exist for all t. Since the defects in our system are planar, the rare regions are 
infinite in the X\ and x-i directions but finite in the remaining directions. This is a 
crucial difference to systems with uncorrelated disorder, where the rare regions are 
finite objects (see Sec. E}. In our system, each rare region is therefore equivalent to a 
two-dimensional Ising system and can undergo a real phase transition independently 
of the rest of the system. Thus, for planar defects, those rare regions which are in the 
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ordered phase will develop a non-zero static order parameter which can be aligned by 
an infinitesimally small interaction or an infinitesimally small field. 

The resulting phase transition is very different from conventional continuous 
phase transitions, be it transitions controlled by clean or finite-disorder or even 
infinite-randomness fixed points. In all of those cases, which represent "sharp" phase 
transitions, a non-zero order parameter develops as a collective effect of the entire 
system, accompanied by a diverging correlation length in all directions. In contrast, 
in a system with planar defects, the order parameter develops very inhomogeneously 
in space with different parts of the system (i.e., different regions) ordering 
independently at different t. Correspondingly, the correlation length in r_L direction 
remains finite for all temperatures. This defines a smeared transition. We thus 
conclude that planar defects lead to a smearing of the phase transition. 

Above, we have defined a smeared continuous phase transition via the behavior 
of correlation length. We now characterize it via the behavior of the free energy 
density or other global thermodynamic variables.! Two qualitatively different cases of 
smeared transitions must be distinguished. If the disorder distribution is unbounded 
in the sense that it permits regions with a local T c — oo (system Q with Gaussian 
disorder falls into this class), the total (average) order parameter is non-zero for all 
temperatures, and the free energy density is analytic. In contrast, if the disorder 
distribution is bounded, there is a true paramagnetic phase with zero order parameter 
at high temperatures. At some temperature To, a non-zero order parameter starts 
to develop on rare spatial regions, accompanied by an essential singularity in the 
free energy density (which stems from the probability for finding a rare region). 
This singularity is weaker than the power-law singularities at the sharp transitions 
mentioned above. The model Q with Poisson disorder falls into this second class, 
with To being identical to the critical temperature of the pure system. 

3. Lifshitz tail arguments 

In this section we use extremal statistics to derive the leading thermodynamic behavior 
in the 'tail' of the smeared transition, i.e., in the parameter region where a few islands 
have developed static order but their density is so small that they can be treated 
as independent. The approach is similar to that of Lifshitz and others for the 
description of the tails in the electronic density of states. We first consider the Poisson 
type of disorder and later describe the differences for Gaussian disorder. 

3.1. Poisson disorder 

The probability w of finding a large region of linear size Lr (in rj_-space) devoid of 
any impurities is, up to pre-exponcntial factors, given by 

w ~exp(-cL^-) . (4) 

As discussed in Sec. |21 such a rare region develops static long-range order at some 
t c (L R ) below the clean critical point t° — 0.§ The value of t c (L R ) varies with the 
size of the region: The largest islands develop long-range order closest to the clean 

\ Note that global thermodynamic variables (which average over the entire sample) do not resolve 
the strong spatial inhomogeneities which are the hallmark of a smeared transition, and thus do not 
provide complete information about the transition 

§ We use the bare value, t® = 0, for the clean critical reduced temperature. For our purpose, the 
distinction between the bare and the renormalized t is of no importance. 
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critical point. Finite size scaling (for the clean system, because the island is devoid of 
impurities) yields 

t° c -t c (L R ) = \t c (L R )\=AL^ (5) 

where <j) is the finite-size scaling shift exponent and A is the amplitude for the crossover 
from d dimensions to a slab geometry infinite in two dimensions but finite in d± = d— 2 
dimensions. Combining and JSJ we obtain the probability for finding a rare region 
which becomes critical at t c as 

w(t c ) ~ exp(-B \t c \-* x l+) (fort -►()-) (6) 
where the constant B is given by B = cA d± ^. The total (or average) order parameter 
m is obtained by integrating over all rare regions which are ordered at t, i.e., all rare 
regions having t c > t. Since the functional dependence on t of the order parameter on 
a given island is of power-law type it does not enter the leading exponentials but only 
the pre-exponential factors. Therefore we obtain to exponential accuracy 

m ~ exp(-B |t c |- dj -/*) (for t -> 0-) . (7) 

The homogeneous magnetic susceptibility consists of two different contributions, one 
from the islands on the verge of ordering and one from the bulk system still deep 
in the disordered phase. The bulk system provides a finite, noncritical background 
susceptibility throughout the tail region of the smeared transition. In order to estimate 
the contribution of the islands consider the onset of local ordering at the clean critical 
temperature t° c = 0. Using JBJ for the density of the islands we can estimate 

,A 

X ~ / dt t^ exp(-B £- f W0) (8) 



which is finite because the exponentially decreasing island density overcomes the 
power-law divergence of the susceptibility of an individual island. Here 7 is the clean 
susceptibility exponent and A is related to a lower cutoff for the island size. Once 
ordered islands exist they produce an effective background magnetic field everywhere 
in space which cuts off any possible divergence. Therefore we conclude that the 
homogeneous magnetic susceptibility does not diverge anywhere in the tail region of 
the smeared transition. However, there is an essential singularity at the clean critical 
point produced by the vanishing density of ordered islands. 

The spatial magnetization distribution in the tail of the smeared transition is very 
inhomogeneous. On the already ordered islands, the local order parameter m(r) is of 
the same order of magnitude as in the clean system. Away from these islands it decays 
exponentially with the distance from the nearest island. The probability distribution 
P[logm(r)] will therefore be very broad, ranging from logm(r) = O(l) on the largest 
islands to logm(r) — > —00 on sites very far away from an ordered island. The typical 
local order parameter mt yp can be estimated from the typical distance of any point to 
the nearest ordered island. From © we obtain 

r tvp ~ exp(B \t\- d ^*/d±) . (9) 

At this distance from an ordered island, the local order parameter has decayed to 

-Cexp(B \t\- d± /+/dx)] (10) 

where £0 is the bulk correlation length (which is finite and changes slowly throughout 
the tail region of the smeared transition) and C is constant. A comparison with (JTJl 
gives the relation between m typ and the thermodynamic order parameter m, 

\logm typ \ ~ m- 1 ^ . (11) 



mtyp ~ e rt yp/i° ~ exp 
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Thus, m t y P decays exponentially with to indicating an extremely broad local order 
parameter distribution. In order to determine the functional form of the local order 
parameter distribution, first consider a situation with just a single ordered island 
at the origin of the coordinate system. For large distances r the magnetization 



falls off exponentially like m(r) = too e r± ^°. The probability distribution of 
x = log[m(r)] = log mo — J"j_/£q can be calculated from 

dN 



log too 

P(\X\): 



dN 


dN 


dr± 


dx 


dr± 


dx 



dr i 



d-l 



(12) 



where dN is the number of sites at a distance from the origin between r± and r± + dr± 
or, equivalently, having a logarithm of the local magnetization between x and x + dx. 
For large distances we have \x\ ~ rx< Therefore, the probability distribution of log to 
generated by a single ordered island takes the form 

P[log(m)] - | log(m)| d - L ~ 1 (for to < 1) . (13) 

In the tail region of the smeared transition the system consists of a few ordered islands 
whose distance is large compared to £o- The probability distribution of log[m(r)] thus 
takes the form i|13|) with a lower cutoff corresponding to the typical island-island 
distance and an upper cutoff corresponding to a distance £o from an ordered island. 



3.2. Finite-size effects 

It is important to distinguish effects of a finite size Lq in the correlated directions 
and a finite size L± in the uncorrelated directions. If L± is finite but Lc is infinite 
static order on the rare regions can still develop. In this case, the sample contains 
only a finite number of islands of a certain size. As long as the number of relevant 
islands is large, finite size-effects are governed by the central limit theorem. However, 
for t — » 0— very large and rare islands are responsible for the order parameter. The 
number N of islands which order at t c behaves like N ~ Lr^ w(t c ). When N becomes 
of order one, strong sample-to-sample fluctuations arise. Using JBJ for w(t c ) we find 
that strong sample to sample fluctuations start at 

/ j \-<P/d± 

l*£j~(%l°g(£0) • (14) 

Thus, finite size effects are suppressed only logarithmically. 

Analogously, one can study the onset of static order in a sample of finite size L± 
(i.e., the ordering temperature of the largest rare region in this sample). For small 
sample size L±, the probability distribution P(t c ) of the sample ordering temperatures 
t c will be broad because some samples do not contain any large islands. With 
increasing sample size the distribution becomes narrower and moves towards the 
clean t® because more samples contain large islands. The maximum t c coincides with 
i° corresponding to a sample without impurities. The lower cutoff corresponds to 
an island size so small that essentially every sample contains at least one of them. 
Consequently, the width of the distribution of critical temperatures in finite-size 
samples is governed by the same relation as the onset of the fluctuations, 

At c ~(2£log(£0) ■ (15) 

We now turn to finite-size effects produced by a finite extension Lc of the system 
in the correlated directions. In this case, the true static order on the rare regions 
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is destroyed. Whether or not the system develops long-range order depends on the 
number d± of uncorrelated dimensions. If d± = 1 no static long-range order can 
develop, i.e. the transition is rounded by conventional finite size effects in addition 
to the disorder-induced smearing. In contrast, for d± > 1 a true phase transition 
is possible but requires a finite interaction between the islands of the order of the 
temperature. This restores a sharp phase transition at a reduced temperature t c (Lc). 
To estimate the relation between Lc and t c (Lc) we note that the interaction between 
two planar rare regions of linear size Lc is proportional to L c and decays exponentially 
with their spatial distance r, Ei nt ~ L c exp(— r/£o), where £o is the bulk correlation 
length. With the typical distance given by © we obtain a double exponential 
dependence between Lc and the critical reduced temperature t c (Lc)' 

log(logL c ) ~ \t B \-*±'+ . (16) 



3.3. Gaussian disorder 

In this subsection we apply Lifshitz tail arguments to the case of 5-correlated Gaussian 
disorder In contrast to the positive definite Poisson disorder, the Gaussian 

disorder is unbounded for 8t(r) —> —oo. Therefore, locally ordered islands can exist for 
all t, and, in principle, the tail region of the smeared transition stretches to t — > oo. In 
order to determine the probability for finding a rare region which orders at a certain t c , 
i.e., the Gaussian equivalent to JBJ, consider an island of linear size Lr. The average 
disorder value (St) in this region has the distribution 

P[(St)]~ exp[-L^(5t) 2 /(2A)] . (17) 

According to finite size scaling, an island of size Lr with an average disorder value 
of (St) will develop static order at t c — —(St) — AL^. For fixed island size, the 
probability distribution of t c reads 

P(t c )~exp[-L d ^(t c + AL-/) 2 /(2A)] . (18) 

We now calculate the stationary point of the exponent with respect to Lr (saddle point 
approximation) to determine which island size gives the dominating contribution. For 
(2 - dx/(f>) > 0, this leads to 

(2-d±/<l>)AL-+=t d/4>. (19) 

This implies that small but very deep fluctuations are responsible for the tail of the 
smeared transition (in contrast to Poisson disorder where the tail is produced by the 
largest islands). Inserting into (|18|) gives the desired probability as a function of t c , 

w(t c ) ~ exp(-B i 2 -^/*) (20) 

where B is a constant. In contrast, for (2 — d±/<p) < 0, the main contribution comes 
from arbitrarily small islands. In this case the ^-correlated disorder is unphysical and 
should be replaced by a distribution with a finite correlation length. The resulting 
w(t c ) is then purely Gaussian. In the following we will not consider this case. 

Starting from (|20[) . one can derive the leading behavior in the tail of the smeared 
transition (t —> oo). Integrating over all ordered islands gives the total magnetization 

m ~ exp(-B . (21) 

The typical local magnetization m t y P decays exponentially with the average 
magnetization m, i.e. lfTT|) also holds in the Gaussian case as can be easily seen by 
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starting the derivation from l|2U|) rather than (jSJ. For samples of finite size in the 
uncorrelated directions the onset of sample-to-sample fluctuations and the width of 
the sample ordering temperatures follow 

t L± ~ At c ~(logL L ) 1/(2 -' W * ) ■ (22) 

If the samples are of finite size in the correlated directions but infinite in the 
uncorrelated ones (and if dj_ > 1) a sharp phase transition is restored. The relation 
between Lq and t c is double-exponential: 

log(logL c )~^ /0 ■ (23) 

A comparison of the results for Gaussian and Poisson disorder shows that the 
functional form of the thermodynamic relations at the smeared transition depends 
on the form of the disorder and is thus not universal. 



4. Numerical results 

4-1- Model and method 

To illustrate the smearing of the phase transition we now show results from a numerical 
simulation of a 3D Ising model with planar bond defects, i.e., dc = 2 and d± = 1. 
As discussed above, the smearing of the phase transition is the result of exponentially 
rare events. Therefore large system sizes are required to observe it. In order to reach 
sufficiently large sizes with a reasonable effort and to retain the possibility of static 
order on an isolated rare region, we consider a model with infinite-range interactions 
in the correlated directions (parallel to the defects) but short-range interactions in the 
uncorrelated direction (perpendicular to the defect planes). While the infinite-range 
model will not be quantitatively comparable to a short-range model, it provides a 
simple example for the rounding mechanism introduced in Sec. |21 Moreover, since 
the rounding mechanism only depends on the existence of static order on the rare 
regions but not on any details, we expect that the results will qualitatively valid for a 
short-range model, too (with the appropriate changes to the exponents in the relations 
derived in Sec. OJ- Indeed, preliminary results |35| of Monte-Carlo simulations of a 
3D short-range Ising model with planar defects fully agree with this expectation. 
The Hamiltonian considered in this section reads 

n jj2 ^ ^ S X y Z S x +i ^' z f ^ ^ Jx^x ,y ,x^x,y' z' ^ ^ ^ ^x,y .x • (24) 

C x,y,z,y',z f C x,y,z,y' ,z' x,y,z 

Here x, y, z are the integer-valued coordinates of the Ising spins. The magnetic 
field h acts as a symmetry breaker, it will be set to a small but non-zero value, 
typically h = 1CP 10 . Lq is the system size in the y and z directions. The planar 
defects are parallel to the (y, z)-plane. They are introduced via disorder in the in- 
plane coupling J x which is a quenched binary random variable with the distribution 
P(J) = (1 — c) 5( J— l)+c S(J— Jd) with Jd < 1. Here, c is the defect concentration and 
the deviation of Jd from 1 measures the strength of the defects. For our calculations 
we use Jd = 0, i.e., very strong defects. The fact that one can independently vary 
concentration and strength of the defects in an easy way is the main advantage of this 
binary disorder distribution. However, it also has unwanted consequences, viz. log- 
periodic oscillations of many observables as functions of the distance from the critical 
point 36 . These oscillations are special to the binary distribution and unrelated to 
the smearing considered here; we will not discuss them further. 
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Because the interaction is infinite-ranged in the correlated (y and z) directions, 
these dimensions can be treated exactly: We introduce the plane magnetizations 



into the partition function and integrate out the original Ising variables S XtV>z . In 
the thermodynamic limit, Lq — > oo, the resulting integral can be solved using saddle 
point approximation. This leads to a set of coupled local mean-field equations 



where T is the temperature (we work in units where the Boltzmann constant is 
ks = 1)- In the clean case, c = 0, all J x are equal to 1. Consequently, all local 
mean-field equations are identical, reading m = tanh[(3m + h)/T]. The clean critical 
temperature of the model H24|) is therefore T c ° = 3. 

In the presence of disorder, the local mean-field equations l|26ll have to be solved 
numerically We employ a simple self-consistency cycle: All m x are initially set to 
m x — 0.1. We then insert the current values of m x into the r.h.s. of (|26[) and calculate 
new values m x cw which are then put into the r.h.s. and so on. This is repeated until 
the change \m x ew — m x \ falls below a threshold of 10~ 12 for all x. This method is 
very robust, it never failed to converge, but it can be very slow in the tail region 
of the transition where ordered islands coexist with large disordered regions. In this 
region sometimes several 1000 iterations were necessary. However, in contrast to 
more sophisticated methods like the Newton-Raphson method no matrix inversions 
are required. Therefore we were able to simulate very large systems up to L± = 10 6 . 

We also note that the model H24|) is similar to the mean-field McCoy- Wu model 
considered in Ref. [37|- The raw data obtained in this paper (L± up to 1000) are very 
similar to ours, but they were interpreted in terms of conventional power-law critical 
behavior with an unusually large order parameter exponent (3 rj 3.6. This may be 
partially due to an unfortunate choice of parameters (in particular, a higher impurity 
concentration of c = 0.5) which makes it harder to observe the exponential tail of the 
smeared transition. We will come back to this question in Sec. El 

4-2. Thermodynamic magnetization and susceptibility 

In this subsection we present numerical results for the total magnetization m = 
LJ 1 J2 X m x an d the homogeneous susceptibility x = dm/dh. The left panel of figure 
^ gives an overview over the behavior of m and x as functions of temperature T for 
size L± — 1000 and impurity concentration c = 0.2. The data are averages over 1000 
disorder realizations. (Thermodynamic quantities involve averaging over the whole 
system. Thus ensemble averages rather than typical values give the correct infinite 
system approximation.) At a first glance these data suggest a sharp phase transition 
close to T = 2.96. However, a closer inspection of the transition region in the right 
panel of figure ^ shows that the singularities are rounded. If this smearing was a 
conventional finite-size effect the magnetization curve should become sharper with 
increasing L± and the susceptibility peak should diverge. This is not the case here, 
both the magnetization and the susceptibility data are essentially size-independent. 
We conclude that the smearing is an intrinsic effect of the infinite system. 

In order to compare the total magnetization data with the analytical Lifshitz 
tail prediction J7J we need to determine the value of the clean finite-size scaling shift 




(25) 



m x = tanh [(m x _i 



+ J x m x + m x+1 + h)/T 



(26) 
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Figure 1. Left: Magnetization m and susceptibility \ f° r Lx = 1000, c = 0.2. 
Right: Magnetization and susceptibility close to the seeming transition for 
different system sizes L± in the uncorrelated direction. 




Figure 2. Left: Total magnetization m as a function of (T c ° -T)" 1 / 2 for several 
impurity concentrations c. The solid lines is are fits to eq. J7J with (j> = 2. Right: 
Decay constant B as a function of — log(l — c). 



exponent for the model l)24|l . Finite-size scaling is governed by the first appearance of a 
non-vanishing solution of the linearized clean mean-field equation in a slab geometry 
of varying thickness. This equation is equivalent to a one-dimensional Schrodinger 
equation in a potential well, leading to a clean shift exponent of 4> = 2 and dj_/4> = 1 /2. 
In the left panel of figure[21 we therefore plot the logarithm of the total magnetization, 
averaged over 300 samples, as a function of (Tj? — T) -1 / 2 for system size Lj_ = 10000 
and several impurity concentrations c. For all c, the data follow eq. (JJJ over several 
orders of magnitude in m. Fits of the data to Q are used to determine the decay 
constants B. The right panel of figure |2] shows that these decay constants depend 
linearly on — log(l — c). This is exactly the expected behavior in a lattice model 
with binary disorder since the probability for finding an island of size Lr devoid of 
impurities behaves like (1 - c) Lr ~ e L « log ( 1 ~ c ). 

4-3. Local magnetization 

In the tail region of the smeared transition the system consists of a few rare ordered 
regions (on sufficiently large islands devoid of impurities). These rare regions are far 
apart, and between them the local magnetization is exponentially small. This behavior 
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Figure 3. Left: Local magnetization m x vs. x for a segment of a single sample 
with size L± = 100000, impurity concentration c = 0.2 and temperature T = 2.99. 
Right: Local magnetization m x and local coupling constant J x in the vicinity of 
one of the ordered islands. 

is illustrated in figure E3 The left panel shows the local magnetization of a segment 
of a large sample in the tail region of the smeared transition. The data show that 
the local magnetization on the islands is large (s=s 0.1) but it drops very rapidly with 
increasing distance from the islands. This drop-off can be used to estimate the bulk 
correlation length which is very small in this parameter region, £o ~ 3... 4. Note that 
the seeming saturation of the magnetization decrease at m « 10~ 9 is a result of the 
finite external field h = 1CP 10 . Without field, the magnetization would drop much 
further. A comparison of the local magnetization m x and the local coupling constant 
J x in the right panel of figure [3] shows that magnetic order only exists on a sufficiently 
large island devoid of any impurities. 

In order to quantify these observations, we have also calculated the probability 
distribution P(logm x ) of the local magnetization values, or, more precise, their 
logarithms. Figure 01 shows this distribution for a single large sample of L± — 200000, 
impurity concentration c = 0.2, and temperatures ranging from T = 2.8 (deep in 
the ordered phase) to T = 2.99 (in the tail of the smeared transition). Clearly, with 
increasing temperatures the distribution becomes very broad, even on this logarithmic 
scale. Two qualitatively different regimes can be distinguished. For temperatures not 
too close to the clean critical point (T = 2.80 to 2.96), the distribution is dominated 
by an exponential decay towards small magnetization values (i.e., towards negative 
logm) with the decay constant decreasing with increasing temperature. However, 
once the system enters the rare-region dominated tail region of the smeared transition 
(which, as can be seen from figures Q an d happens roughly at T = 2. 96. ..2. 97 for 
c = 0.2) the behavior changes: In agreement with (|13|l the probability distribution 
P(logm x ) becomes constant except for a peak at large magnetization values which we 
attribute to the behavior of m x on the ordered islands (eq. (|13fl was derived from the 
behavior of the sites not belonging to one of the rare ordered islands). 

4-4- Finite-size effects and sample-to-sample fluctuations 

In this subsection we present numerical data for the effects of a finite sample size L±_ in 
the uncorrelated direction, i.e., perpendicular to the impurities. Effects of a finite size 
he in the correlated directions cannot easily be studied within our approach because 
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Figure 4. Probability distribution of the local magnetization m x for single 
sample of size L± = 200000, impurity concentration c = 0.2 and several 
temperatures. The three curves represented by filled symbols (T > 2.98) are 
in the asymptotic regime described by eq. 1131 . 
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Figure 5. Left: Average and typical (log. averaged) sample magnetization 
for c = 0.2 and different sample sizes. The inset shows the onset temperature of 
strong sample-to-sample fluctuations (measured by = — T^) as a function of 
(logL^) - 2 . Right: Probability distribution of the sample ordering temperatures 
for c = 0.2 and different sample sizes. The inset shows the width AT C of this 
distribution as a function of (logLx) -2 - 



the limit Lq — > oo has been taken in deriving the local mean-field equations 126|) . 

To investigate sample-to-sample fluctuations in a finite-size system we compare 
the average sample magnetization m av = (m) and the typical sample magnetization 
TTityp = exp(logm) where m is the total magnetization of a sample and (■) denotes 
the average over the disorder. A significant difference between average and typical 
magnetizations indicates strong sample-to-sample fluctuations. The left panel of figure 
[SJshows the typical magnetizations for four different sizes L± together with the average 
magnetization (which is size-independent within the accuracy of this plot) for 300 to 
1000 samples, depending on L±. In the tail region of the smeared transition the typical 
and average magnetizations deviate from each other by several orders of magnitude. 
From these data we define a temperature Tl ± which describes the onset of the sample- 
to-sample fluctuations via m av (TL ± )/mt yp (TL ± ) = e. This temperature shifts to larger 
values with increasing system size L±. The inset shows tL ± = T" c ° — Tl ± as a function 
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of (logLjJ~ 2 . The plot shows that tj JX follows eq. 11411 in good approximation. 

The right panel of figure[5]shows the probability distribution of the sample critical 
temperature which is defined as the ordering temperature of the largest rare region in 
a particular sample (600 to 2000 samples used). As predicted in subsection 13.21 the 
distribution becomes narrower with increasing size L± , and it shifts towards the clean 
critical temperature T c ° = 3. The width AT C of this distribution has been determined 
from the values where P has decreased by a factor of 1/e. The inset shows that AT C 
is proportional to (logi^)~ 2 as predicted in eq. I)15[l. 

5. Discussion 

To summarize, we have shown that true static order can develop on an isolated rare 
region in an Ising model with planar defects. As a result, different parts of the system 
undergo the phase transition at different temperatures, i.e., the sharp transition is 
smeared by the defects. In this final section we discuss the generality of our findings 
and their relation to Griffiths phenomena and the Harris criterion. We also discuss 
favorable conditions for observing the smearing in experiments or simulations. 

The origins of the smearing of the phase transition introduced in this paper and 
of Griffiths phenomena are very similar, both are caused by rare large spatial regions 
which are locally in the ordered phase. The difference between the two effects is a 
result of disorder correlations. For uncorrelated or short-range correlated disorder, the 
rare regions are of finite size. Thus, they cannot develop true static long-range order. 
Instead, the order parameter fluctuates slowly leading to the Griffiths singularities [H] 
discussed in Sec. ^ In contrast, if the rare regions are infinite in at least two 
dimensions, a stronger effect occurs: the rare regions which are locally in the ordered 
phase actually develop true static order, and this leads to a smeared transition. In 
other words, the same rare regions which would usually produce Griffiths singularities 
are responsible for the smearing of the transition in a system with planar defects. 
The tail of the smeared transition stretches to where the Griffiths temperature would 
normally be, i.e., to T — oo for Gaussian disorder or T — Tj? for Poisson disorder. So, 
the smeared transition replaces not only the conventional critical point but also the 
Griffiths phase. 

We now discuss under which conditions a phase transition will be smeared by 
disorder. The arguments given in Sec. |21 in favor of a smeared transition only relied 
on the dimensionality of the rare regions and the order parameter symmetry. It is 
therefore expected that all phase transitions with discrete order parameter symmetry 
are smeared by quenched disorder if the defects are perfectly correlated in at least 
two dimensions. In systems with linear defects the transition will remain sharp 
if the interactions are short-ranged. However, if the interactions in the correlated 
direction fall of as 1/r 2 or slower, even linear defects can lead to smearing because a 
ID Ising model with 1/r 2 interaction has an ordered phase [SHj. This is particularly 
important |39j for quantum phase transitions in itinerant electronic systems which can 
be mapped onto classical systems with 1/r 2 interaction in imaginary time direction. 

Phase transitions with continuous order parameter symmetry are more stable 
against smearing. In systems with short-range interactions the dimensionality of the 
defects has to be at least three, and in systems with linear or planar defects long- 
range interactions are required for smearing. It is known (40j that classical XY and 
Heisenberg systems in dimensions d = 1, 2 develop long-range order at finite T only if 
the interaction falls off more slowly than l/r 2d . Consequently, in a system with linear 
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defects the phase transition will only be smeared if the interactions in defect direction 
fall off more slowly than 1 /r 2 (or more slowly than 1 /r 4 for planar defects) . 

A third important question concerns the universality of the thermodynamic 
behavior in the vicinity of the smeared transition. As can be seen from comparing 
the results of Poisson and Gaussian disorder, the functional dependence of the 
magnetization and other observables on the temperature is not universal, it depends on 
details of the disorder distribution. Therefore, only the presence or absence of smearing 
is universal in the sense of critical phenomena (i.e., depending on dimensionality and 
symmetry only) while the thermodynamic relations are non-universal. 

Our fourth remark deals with the relation of the disorder-induced smearing and 
the Harris criterion which in the case of extended defects reads [21] v > 2/d±. We 
emphasize that the phase transition can be smeared by planar defects even if the 
corresponding clean critical point fulfills the Harris criterion and appears to be stable. 
The reason is that the Harris criterion assumes a homogeneous transition and studies 
the behavior of the coarse-grained (root-mean-square) disorder at large length scales. 
However, the formation of static order on an isolated finite-size rare region is a non- 
perturbative finite-length scale effect in the tail of the disorder distribution. This type 
of effects is not covered by the Harris criterion. 

We now turn to the question which conditions are favorable for the observation of 
the smearing in experiments or simulations. A significant number of large rare regions 
will only exist if the defect concentration is small, as can be seen from (JBJ. On the 
other hand, the impurities have to be sufficiently strong so that the bulk system is 
still far away from criticality when the first rare regions start to order. Thus, the most 
favorable type of disorder for the observation of the smearing is a small density of 
strong impurities. In contrast, for weak impurities the rounding is restricted to a very 
narrow temperature range below the clean critical temperature and maybe masked 
by the remainder of the bulk critical fluctuations. If the impurity concentration is 
too high, the exponential drop-off is very fast making it very hard to observe the 
smearing. This may also explain why no smearing was observed in earlier simulations. 
Specifically, in Ref. [23] weak impurities (A J/ J = 0.1) of a high concentration c = 0.5 
were used. These are unfavorable conditions, and the maximum system size of L = 27 
used in this simulation was probably too small to observe the smearing. 

Finally, we add one remark about the extreme tail of the smeared transition. In 
this parameter region, the ordered islands are very far apart and the aligning effect 
of the direct island-island interactions is very weak. It is therefore possible that in 
a real system other interactions play an important role in determining the relative 
orientation of the order parameter on the islands. This could lead, e.g., to spin glass 
order in the extreme tail of the smeared transition. 

In conclusion, perfect disorder correlations in one or more directions dramatically 
increase the effects of the disorder on a phase transition with discrete order parameter 
symmetry. Zero-dimensional defects typically lead to a conventional critical point with 
power-law scaling. For linear defects the generic behavior seems to be an infinite- 
randomness critical point 24 , 25 2£f| while we have shown here that planar defects 
generically destroy the sharp transition by smearing. 
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